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Abstract 



Recently D. Vollick [Phys. Rev. D68, 063510 (2003)] has shown that the inclusion of 
the 1/R curvature terms in the gravitational action and the use of the Palatini formalism 
^ , offer an alternative explanation for cosmological acceleration. In this work we show not 

0^ I only that this model of Vollick does not have a good Newtonian limit, but also that any 

^ ■ f{R) theory with a pole of order n in R = and {d? f /(PR){Rq) 7^ 0, where Rq is the 

, scalar curvature of background, does not have a good Newtonian limit. 

o 

g . PACS number(s): 95.35.+d, 04.25.Nx, 98.80.-k, 04.20.-q 

"o 

7: 1 INTRODUCTION 



From recent studies it seems well established that our universe expansion is currently in an 
accelerating phase. The evidence of cosmic acceleration has arisen not only from the high 
redshift surveys of type la supernovae [1, 2], but also from the anisotropy power spectrum of 
the cosmic microwave background [3, 4]. One of the most accepted explanations is that the 
Universe has been dominated by some form of dark energy for a long time. However, none of 
the existing dark energy models are completely satisfactory. 

It is possible to find other explanations for cosmic expansion using field equations other than 
Eisntein's equations. Recently, some authors have proposed to add a R'^ term in the Einstein- 
Hilbert action to modify general relativity [5, 6, 7]. They obtained the field equations using 
second-order formalism, varying only the metric field, and thus obtained the so-called fourth- 
order field equations. Although the models were obtained using corrections of the Einstein- 
Hilbert Lagrangian of type i?", where n can take a positive or negative value to explain both 
the inflation at an early time and the expansion at the present time [5, 8], they still suffer from 
violent instabilities [9]. The Newtonian limit of these fourth- order theories has been studied 
by Dick [10]. 

* Electronic addresss: domingue@fis.uncor.edu 
'^Electronic address: barraco@famaf.unc.edu. ar 



1 



On the other hand, we can consider those theories that are obtained from a Lagrangian den- 
sity Ct{R) = f {R)y/--g + Cm, which depends on the scalar curvature and a matter Lagrangian 
that does not depend on the connection, and then we can apply Palatini's method to obtain the 
field equations [12, 13]. In Refs. [12, 13], we showed the universality of the Einstein equation 
using a cosmological constant. More recently, Ferraris, Prancaviglia and Volovich pubhshed 
the same result [14, 15]. For these theories we have studied the conserved quantities [16], the 
spherically symmetric solutions [17], the Newtonian limit [12, 18], and the cosmology described 
by Friedmann- Robertson- Walker (FRW) metric [19]. 

In our previous paper [18], it was also shown that it is very difficult to test these models in 
the (post-) Newtonian approximation. The reason for this is that the departures from Newtonian 
behavior arc both very small and are masked by other effects, because these departures from the 
Newtonian behavior have to be measured when the body is moving "through" a matter-filled 
region. In the main applications of these theories, for example, in the Newtonian limit or in 
the calculation of the observed cosmological parameter, we have considered the analyticity of 
f{R) at = 0. 

Moreover, it is well known [18, 14, 17] that in a vacuum, or in the case of T = const. 
the solutions of these theories are the same as general relativity with a cosmological constant, 
even when f{R) is not analytical at i? = 0. On the other hand, solutions corresponding to 
different cosmological constants are allowed by some of these theories. Therefore, the homoge- 
neous and isotropic vacuum solution for these theories is the de Sitter space-time with different 
cosmological constants, except when one of the allowed cosmological constants is A = 0, which 
corresponds to flat space-time. 

Recently, Vollick [20] used the corrected Lagrangian of the works of Carroll et al, and 
Capozziello et al. [5, 7], f{R) = R — a'^/R, and the Palatini variational principle -which is a 
particular case of the above theories but is not analytical at i? = 0- to prove that the solutions 
of the field equations approach de Sitter universe at a late time. This result was obtained 
using the above well known property of the vacuum solutions. Thus, the inclusion of l/i? 
curvature terms in the gravitational Lagrangean provides us with an alternative explanation 
for the cosmological acceleration. Moreover, Meng and Wang [21, 22] have studied the modified 
Friedmann equation with the Palatini variational principle and its first-, second-, and third- 
order approximated equations. 

In agreement with this last idea about the maximal symmetric vacuum solution, we have 
studied the Newtonian limit for f{R) theories as a perturbation of the de Sitter background. 
The difference from previous works [12, 18] is that there the background was fiat space-time, 
which is only possible in the particular case /(O) = 0. 

In this work we show that the Vollick model [20] does not give the Newtonian limit but 
also that any f{R) with a pole of order n in R = and {(P f / (P'R){Rq) ^ 0, where Rq is the 
scalar curvature of background, cannot have a Newtonian Limit. This result is equivalent to 
the condition obtained by Dick [10] for the fourth-order theories. Nevertheless, the theories 
with singular f{R) and [d? f / cP R){Rq) — must be studied in the future since they have a 
correct Newtonian limit and explain the observed cosmological acceleration. In our work we 
follow the conventions and notation of Synge [11] 
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2 REVIEW OF THE FIELD EQUATIONS 



We review in this section the structure of the theory as presented elsewhere in previous works 
[12, 18]. 

Let Ai he a manifold with metric gab and a torsion-free derivative operator Va, both con- 
sidered as independent variables. Consider a Lagrangian density C — f{R)y/—g + jOm, where 
the matter Lagrangian Cm does not depend on the connection. 

Suppose we have a smooth one-parameter A family of field configurations starting from 
given fields g""^, Vo, and t/j (the matter fields), with appropriate boundary conditions. Let 5g"'^, 
Sr^f^, Sip be the corresponding variations of those fields, i.e., Sg"-^ — (dg'^/dX) \x=o, etc. Then, 
if we vary with respect to the metric, the field equations are 

f'{R)Rab - \f{R)gab = OiMTab. (1) 

where f'{R) — (df/dR), {SSu/Sg"'^) = —Tab^/—g and = — Stt. The variation with respect 
to the connection, recalling that this is fixed at the boundary, gives 

VA^gg'-'fm = 0. (2) 

Now, we choose Lagrangian f{R) with f'{R) derivable. Then, the last equation becomes 

^cgab = Wgab, (3) 

where 

6e--[ln/'(i?)],e. (4) 

Thus, we have a Weyl conformal geometry with a Weyl field given by Eq. (4). 
Prom Eq. (1) we obtain 

f{R)R-2f{R)^aMT, (5) 

which defines R{T), and we suppose that the function f{R) is such that R{T) is derivable with 
respect to the variable T. Therefore, is determined by T and its derivatives except in the 
case f{R) = ujR^, for which Rf — 2/ = 0, so we must consistently have T = 0. It is important 
to note that b^ has a solution only if T is differentiable in 7W; this condition on T, for the 
existence of a solution, is not necessary in other theories such as general relativity (GR) or 
fourth-order theories. 

The connection solution to Eq. (3) is 

r^c = Cl -\{5tb, + S^^b, - g,r) , (6) 

where C"" he are the Christoffel symbols (metric connection). Then, we have to solve only Eq. 
(!)• 

The Riemann tensor can be defined in the usual form, and then, the Ricci tensor and scalar 
curvature are 

Rab = Rab - \Dabh + ^D^ba - \gabD ■ b - Uah + ^^ab&' (7) 

R^R^ -?,D-b+h)\ (8) 
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where i?^, R^, and Dc are the Ricci tensor, scalar curvature, and covariant derivative defined 
from the metric connection, respectively. 

Because the matter action must be invariant under difi^eomorphisms and the matter fields 
satisfy the matter field equations. Tab is conserved [18] 

D-Tab = 0. (9) 

Therefore, a test particle will follow the geodesies of the metric connection. Using Eqs. (4) and 
(5) we have 

Except for the case of GR, where /" = 0, the Weyl field is nonzero wherever the trace of the 
energy-momentum tensor varies with respect to the coordinates. If T is constant, then R is 
also constant, be = and (5) takes the form 

Gab - l^^Qab = KTab, (11) 

where A and K are two functions of R. All those cases with constant trace of the energy- 
momentum tensor are equivalent to GR for a given cosmo logical constant. This is the so-callcd 
universality of the Einstein equations for matter for which T is constant [12, 13]. In the case 
T — 0, the scalar R is any of the roots, Ri, of the equation f'(R)R — 2f{R) — 0. For each 
root the solutions of the field equations are the solutions of GR with cosmological constant 
A = —Ri/A. Therefore, the maximal symmetric vacuum solution of these theories is the de 
Sitter space-time. 



3 THE NEWTONIAN LIMIT 

The above theories, which explain the cosmological acceleration, must be checked in the 
Newtonian limit, i.e., in the slow field and slow velocity approximation. In this approximation, 
these theories must be in agreement with the present data. 

The Newtonian limit, in our case, must be taken as a perturbation of the homogeneous 
isotropic vacuum background. As we have seen above, this background space-time is the de 
Sitter space-time, except when A = 0, in which case we must perturbate around the fiat space- 
time. The case A = has just been studied [12, 18] and the general solution of the problem 
with A 7^ must agree with the previous result when A = 0. 

The background metric and the background Ricci tensor are 

9ab= -de + e^'^\dx^ + dy^ + dz^), (12) 

Rab^ -A Qab . (13) 

We assume that the metric can be written in the form 

gab ^9ab +hab (14) 
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where hab is the perturbation of the metric. The first order of the field equations (1) are 

/' (R) Rab + /' {R) Rab -\ f (R) 9ab -\ f {R)hab = «m Tab (15) 

where the superscript means zero-order quantities, i.e., background quantities, and the su- 
perscript 1 means the first-order quantities. We have to calculate each term of the above 
equation. 

According to (6) the connection F^^ can be spht into two parts, 

r^, J^t, + hb, (16) 

1 

where V^ab is the connection corresponding to the background metric and V^ab has two terms, 
one depending on the metric perturbation hab a-nd the other depending on which has an 
order higher than zero. 

In order to be able to write the linear field equations we must calculate the linear Ricci 
tensor by using (16). We obtain 

Rab ^Rab + Rab (17) 

1 

where theRab part is 

Rab^ - + V(ar%)c ■ (18) 



By V we mean the covariant derivative associated with the background metric. 
As we have just said, the first order of connection can be written as follows: 

T\b = C'^ab + i^fe, (19) 



where the first term is 



and the second term is 



1 O n I O n 

C^b^g'^Via hb)d - 2 5 hba (20) 

A'ab^ -^iSl lb +51 ha - 9abh- (21) 



By substituting (19) in (18) we obtain 



1 1^ 1^ 1^ 1^ 

Rab='V{aC'^b)c — VcCab + V(a^''fe)c " Vc^'^ab ■ (22) 



In turn, straightforward calculations lead us to 



1 1 ^ ^ 

V(bC''a)c= 2 VaVfe h, (23) 

1 In 

Vc C"afe=V(aV" hb)c+ Rlba) hec' Re{a hb)e hab- (24) 
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Hence, by replacing Eqs. (23) and (24) in (22) we obtain 

1 loo 00^ 0^ ^ 0^ 

Rab= 2 VaVb h— V(aV'^ hb)c— Rfbaf ^ec+ -R''(a ^6)e 



4 ft. 



01 1,01 01 OOl 
ab - 2 V{a66) +^{Vabb + Vbba - dabVch"). 



(25) 



As is well known, there is a gauge freedom in any geometrical theory of gravitation cor- 
responding to the group of diffeomorphisms of space-time. In practice, these diffeomorphisms 
may be viewed as coordinate freedom which may be used to impose coordinate conditions. For 
instance, we may choose coordinates so that, in the linear approximation, the perturbation 
hab and the vector field ha satisfy the gauge condition 



1 1 

he --Vbh+ bb= 0. 
In this gauge the linearized Ricci tensor simplifies to become 



1 In 0^ ^ 0^ 1 1^ 

Rab= 7^ ^ hab— Rlba) ^ec+ -R^(a ^fe)e " 7^ 9abVjf 



(26) 



(27) 



2 ■-u,<j -''(6a) ' -- ^^c^ -uyc 2 

Finally, to rewrite the hnear field equations (15) we have to take into account the factors 

1 

/' {R) and / {R) which depend on T. From (17) we can split the curvature scalar into a zero 
order part and a first order part: 

R = i?o+ k (28) 
where Rq — — 4A. Therefore, using Eq. (5) for first terms we can obtain i? as a function of T: 



1 

R- 



C^M T 



f"{Ro)Ro - f'{Ro) 



(29) 



From this result we can easily obtain / (R) and /' (R). 

Now we are ready to rewrite the first order field equations, in the gauge (26), using Eqs. 
(10), (15), (27), (29); and (13): 



I hab 







4A A 1 

~^hab + -^h 9ab -- 



aM f" 



n 1 
— 9ab^T 

f if" 4A+ /') 







I 



2(4 f A+ h 



1 + 



v 





2r A 





1 

T9ab 



(30) 



It is important to note that when A = we recover the result obtained in Refs. [12, 18]. 

In the Newtonian hmit, the equation of motion of a test particle is given by (9) with v!^ ~ , 
and the proper time of the particle may be approximated by the coordinate time, t. Since the 
sources are " slowly varying" , we expect the space- time geometry to change slowly as well, so 
that the time derivative can be assumed to be negligible. Thus, we find 



00,/i, 



/. = 1,2,3. 



(31) 



As a consequence, we have to study the Newtonian hmit of /iqo using Eq. (30). Taking into 

account the expressions for fl /iqo and At, and the gauge condition (26) in order to substitute 
^12,2, we obtain the field equation for hoo 



g-2t^Ap y2 



Ap 



B 



where 



OiM f" 



■ 

2 /' (4 /" A+ /') 







l\2\ 



C 



(8 A f"f'+if 

0' 

(f)2(4A /" + /') 



(32) 



(33) 



(34) 



Finally, the candidate to be the Newtonian potential is 



■i\x 



X 



X — X' 



d-^x' + Ap{x). 



(35) 
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According to the observed cosmological acceleration, we have to choose A ~ 10" 
Therefore, the last result is consistent with our previous assumption about the time derivative. 
Moreover, the first term in (35) behaves as the Newtonian potential for any current experiment 
or observation. 

In order to ignore the second term, we have to assume | A/C \<^ 1. After a little algebra, 
this inequality can be rewritten as 



r{Ro)/f{Ro) |« 



1 



l + 2Ro 



(36) 



When f{R) is analytical in i? = we can fulfill this condition without problems since Rq ~ 0; 
but we are interested in the singular case, i.e. ,when f{R) is singular in i? = 0. 

If f{R) has a pole of order n in i? = 0, the leading order singularity of f"{R) and f'{R) are 
i?""^^ and respectively. Therefore, the inequality (36) cannot be fulfilled except for the 

particular case 

r (i?o) = 0. (37) 

The inconsistency of the Newtonian limit cannot be attributed to the gauge selection since 
the former is already present in the perturbation equation (15), which was derived previous to 
choosing the gauge. In reality, the problem is originated by ha whose first order of perturbation 
around _Ro has the factor f"{Ro)/f'{Ro). In the singular case this factor becomes greater as Rq 
gets smaller. 

It is interesting to notice that the condition /"(-Rq) = is the same as Dick's condition 
[10] for the case of the singular fourth order theories considered by him. This similarity is not 

surprising on account of the fact, which we will prove elsewhere, that if f{R) = R + ag{R), then 
for the first-order of a the field equations of the fourth order theories and the field equations 
of the Palatini theories are the same. 
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It is not difficult to show that if the Lagrangian density R — a/3R does not satisfies the 
condition (37), then the perturbation around the maximally symmetric vacuum solution does 
not give a Newtonian limit. Nevertheless, as Dick has shown, there exist Lagrangian densities 
that satisfy this condition, for example, f{R) = R + i?^/9/i^ — 3n'^/R. 

4 CONCLUSIONS 

We have studied the first order of perturbation from the maximally symmetric vacuum solutions, 
de Sitter space-time, in a family of theories which are obtained by using the Palatini formalism 
on a general Lagrangian f{R). 

We have proved that the first order of perturbation does not give a good Newtonian limit 
when f{R) has a pole of order "n" in R — 0, and f"{Ro) 7^ 0. However, except for a small 
correction proportional to p, the perturbation in the analytical case has a good Newtonian 
limit. In this analytical case we have also recovered the result of Refs. [12, 18] when A = 0. 

While we were working on this paper Meng and Wang [23] claimed to have shown that the 
Newtonian limit is obtained not only for analytical Lagrangian densities but also for the singular 
case. However, we disagree with their results. We suppose that the main problem is that they 
were not consistent with the first-order of perturbation, since they introduced terms VbVaf'{R) 
instead of using the first order forms of these terms, namely, /"(i?o) / f'{Ro)[Rof"{Ro)^f'{Ro)]^b^aT- 
Another point of disagreement, but not essential to show that the Newtonian limit has prob- 
lems, is about the gauge. They chose the traceless gauge but this gauge does not exist when 
there are sources. 

Finally, the Lagrangian of CaroU et al, and Capozziello et al. [5, 7] does not satisfy the 
condition f"{Ro) = 0, and thus we have proved that the theory of Vollick [20] is not in accor- 
dance with the experimental data. Nevertheless, it is worthwhile to study the f{R) theories in 
the Palatini formulation with an f{R) singular when they fulfill condition (37), because they 
can give a good explanation of the observed cosmological acceleration. 
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